Edge Transport in 2D Cold Atom Optical Lattices 
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We theoretically study the observable response of edge currents in two dimensional cold atom 
optical lattices. As an example we use Gutzwiller mean-field theory to relate persistent edge currents 
surrounding a Mott insulator in a slowly rotating trapped Bose-Hubbard system to time of flight 
measurements. We briefly discuss an application, the detection of Chern number using edge currents 
of a topologically ordered optical lattice insulator. 
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The response theory of transport is a remarkably pre- 
cise framework used to analyze the observable effects of 
applied potentials in a broad class of solid state systems. 
It is natural to ask how experiments on neutral cold 
atoms confined to optical lattices!!]] , predicted to hold a 
variety of novel phases of matter [2| similar to those found 
in the solid state, can make contact with an equivalent 
quantitative framework. Recent experimental work on 
cold atom optical lattices demonstrates essential ingre- 
dients in establishing quantitative response: applied po- 
tentials and detection of conserved quantities. A primary 
tool for detection relies on time of flight (TOF) imaging 
which provided the first evidence for Bosc-Einstein con- 
densation [3j and revealed phases of optical lattice real- 
izations of Bose-Hubbard (BH) models 0, H, @, 0, S] • 

By combining TOF with externally applied potentials 
recent work has demonstrated transport in one dimen- 
sional optical lattices [3, E3] • In a closed two dimensional 
system the notion of "transport" is less direct. A recent 
experiment has applied rotation to weak lattices ll| con- 
fining bosons. While far from the single band BH limit, 
this experiment reveals vortex pinning arising from the 
weak lattice. Recent work [12, El also suggests that 
uniform effective magnetic fields (equivalent to rotation) 
may be applied to optical lattices already in the BH limit. 
Either implementation, rotation or an effective magnetic 
field, can be used as an applied potential valuable in es- 
tablishing persistent currents, and therefore transport, in 
two dimensional lattices. 

Concurrent with experimental progress, a variety of 
cold atom phases have been proposed in two dimensional 
optical lattices 0. Some of the proposed lattice mod- 
els have rich phase diagrams with particularly intriguing 
or even unkno wn g round states, including: extended BH 
models @, El EI, higher band spin models frac- 
tional quantum Hall models [HI], and the Kitaev spin 
model [13, EH EH . We ask how insulating phases arising 
in two dimensional lattice models can be studied using a 
combination of externally applied potentials and TOF. 

Below we argue that trapping leads to edge states 
which serve as a probe of bulk insulating states. As a 
concrete and relevant example we study the slowly rotat- 



ing BH model in detail. Other studies have considered 
vortex configurations in the superfluid phase of the ro- 
tating uniform BH model 2(3, 21 1. Here we study edge 
effects in the Mott insulating phase of the slowly rotat- 
ing trapped BH model. We propose that diamagnetic 
response of edge states can indeed be observed thereby 
offering a quantitative response probe of a variety of bulk 
two dimensional insulators. We briefly discuss implica- 
tions for another insulator where edge states can be used 
to detect the Chern number [It], [HI, [23| of a topologi- 
cally ordered insulator, the non-Abelian ground state of 
the Kitaev model. 

We first note that response to externally applied fields 
can be obtained at a quantitative level by analyzing 
TOF measurements. TOF can be related to the mo- 
mentum density, pk, of particles with lattice momen- 
tum k originally trapped in an optical lattice. Obser- 
vation of pk (with sufficient accuracy) can be combined 
with input parameters to restore quantities of the form: 
J = J2kWkPk, where Wfe is any function of fc which 
can be accurately determined from input experimental 
parameters. By defining Wk = Mk{dEk/dk a ) we ob- 
tain two examples: the free-particle number and energy 
currents in the direction a with the choices Aik — 1 and 
■Mk = Ek, respectively. Here, Ek is the single particle 
energy determined by optical lattice parameters. As we 
will show, diamagnetic current flowing along sufficiently 
narrow edges of optical lattice insulators can be written 
in the form J allowing restoration of the edge current. 

To study the observable response of insulating states in 
trapped optical lattices we consider the two-dimensional 
BH model on a square lattice in the presence of rotation 
(or, equivalently, an effective magnetic field 13, IH) as a 
first step in establishing quantitative response in systems 
nearest ongoing experiments. Using the Peierls substitu- 
tion [3] the BH model in the rotating frame is: 



H = -t jexp {iAi > j)a\a i + h.. 
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where a\ and are the boson creation and number oper- 
ators at the site i, respectively. The parameters include 
the hopping, t, the onsite interaction energy, U, and the 
chemical potential, p. The last term is due to the trap- 
ping potential which adds a site dependent chemical po- 
tential at the square lattice coordinate in the xy plane, 
T{ = (i x ,i y ), parameterized by a modified trapping pa- 
rameter k = k — m(Oa) 2 /2, in units of the lattice spacing, 
a (half the wavelength of the lasers defining the optical 
lattice). The trapping parameter is modified by a term 
due to rotation, with angular frequency, O, of particles 
of mass to. In what follows we find that, for Rb 87 atoms 
with n/U = 0.4, U/E R ~ 0.1, and k/U = 1.2 x 10~ 3 
the modification is negligible for the rotation frequencies 
studied here giving k ks k. En = h 2 /(8ma 2 ) is the recoil 
energy. The rotation also modifies the hopping term to 
give a phase: Ai,j = (■k 2 MI/2Er) J? (z x r) -ds, which can 
be thought of as an integral over a vector potential due 
to an effective magnetic field, B*z, acting on an effective 
charge q* such that q*B* = 2rMl. In the rotating frame 
the neutral bosons experience an effective magnetic field 
which induces a persistent current opposing the applied 
effective field. 

In the linear response regime, Ml /En C 1, we can 
define [25| a number current in the direction a in 
response to the static applied potential, A: Ji^ a = 

ithr 1 ^a|a i — h.c.^j — h~ l Ai^Ki^ where the first term 
is the paramagnetic number current and the second term 
is the diamagnetic number current which contains the ki- 
netic term: Kij a = t ( a\a i + h.c.^ . The total diamag- 
netic number current is an observable response to our ex- 
ternally applied field giving: J® = —h~ 1 '^2 ll Ki^ a Ai : i a . 
In general J D cannot be written in the form J and it 
is therefore not clear how we can relate such a quan- 
tity to TOF measurements. In what follows we will 
show that when the diamagnetic current is confined to 
the edge of the system (and therefore A varies suffi- 
ciently slowly) we can write an approximation to J D 
which in turn can be related to TOF. Consider the follow- 
ing approximation: = —h" 1 Ai^^i J^i-^Ma) where: 
= {h 2 MI/2Er){z X Tie) ■ 6. Here r,e indicates the 
average position of the edge superfluid order parameter, 
xjji = {a}), giving r ■ n* = {J2i \A\ 2 f ■ n)/(E l \fi?)- By 
Fourier transforming and taking the expectation value 
with respect to the ground state we find: {Jj/*) = 
— 2tAi* > ie (HZ) -1 s Pk cos (fc • ft) i where Z is the num- 
ber of nearest neighbors with lattice vectors 8. We now 
have a quantity written in terms of the lattice momentum 
distribution: p k = £V . exp(— ifc • (r, - *v))( i°j)> which, 
we assert, yields an accurate measure of the diamagnetic 
current provided the current flows along the edge. Our 
assertion can be written: (J<P) e « (>^p)e, where () e indi- 
cates averaging in a ground state with only edge current. 
As we will see this relation allows us to probe the edge 
flow around bulk insulators in optical lattices but does 



not necessarily hold for the rotation of a bulk superfluid 
in a trap. To continue with our example of the rotating 
BH model we relate {Jg*) to an observable TOF signal. 

TOF signal can be directly related to the momentum 
distribution even in a slowly rotating optical lattice. 
In the following we assume that the particles do not 
interact after release from the trap. We may then 
apply the free particle propagator K p to a single 
particle Bloch state in the rotating frame, <fik( r )i 
initially confined to the lattice. We project it onto 
a imaged state $ s with imaged coordinates r s in the 
laboratory frame. For slow rotation we find: $ s (r s ) = 
J K p (r;r s )(j)k(r)dr oc S' k q\w (k — AtQz x fe) | 2 , where 
At is the time taken to propagate from the lattice to 
the imaging screen, §' indicates equivalence up to a 
reciprocal lattice vector, and w is the Fourier transform 
of the non-rotating Wannier function. Here the lattice 
wavevector gets mapped to position on the screen in 
free particle propagation: Q{r s ) = (mr s )/(hAt). We 
have derived the above expression to lowest order in 
(hil/Eji) 2 , consistent with our linear response approx- 
imation. The imaged total density is then: n s (r s ) w 
(m/(hAt) f p Q(rs) \w(Q{r s ) - Attl'z x Q(rJ)\ 2 . We 
have found, as in the non-rotating case [26|, [27[, that 
up to an overall Gaussian- like function, w), the imaged 
density on the screen gives p k - We now study the slowly 
rotating BH model under the assumption that pk can be 
accurately extracted from measurements. 

We calculate the ground state of the rotating BH model 
using a modification of the Gutzwiller mean- field ansatz 
0, Ha| . We assume a product state in the Fock number 
basis \Ni), of the form: * = fJi J2n°=o fm \ N i)' where 
the (N c + 1)N complex variational parameters p N . are 
chosen to minimize the ground state energy of H on N 
lattice sites. In what follows we choose N = 50 x 50 where 
the confinement forces the atoms to occupy no more than 
w 45 x 45 sites. We also find that N c = 5 gives suitable 
convergence for the low chemical potentials studied here. 
We minimize (H) using the conjugate gradient method. 
To treat large systems we have developed a three step 
minimization procedure with a computational cost that 
scales linearly with N. Using our product ansatz we first 
find the ground state assuming that each site is an in- 
dependent system with O = 0. In our second step we 
minimize the energy of the whole system using step one 
as an initial guess, while keeping O = 0. This step shows 
29 1 excellent agreement with Monte Carlo simulations 
27| . In the third step we take the variational parame- 
ters of the non-rotating system and modify them to gen- 
erate an initial guess for the rotating system. We use: 

/Arjinitiai = exp (i9iNiV) (l+r?Vi)/wJ n=0 > wnere the ad- 
ditional variational parameter, V , is an integer, rj is a 
random complex number, and 0i is the angular coordi- 
nate of the site i. The above ansatz introduces a vor- 
ticity, V , while finite r\ ensures that our minimization 
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FIG. 1: Total static kinetic energy of the rotating trapped 
Bose-Hubbard model in the rotating frame plotted versus ro- 
tation frequency for parameters (t/U = 0.03 and fi/U = 0.4) 
giving a bulk Mott insulator surrounded by edge superfluid. 
The steps correspond to increasing vorticity V = — 3. The 
top panels plot the momentum distribution in the k x k y plane 
for hfl/E R = 3.31 x 10" 3 (left) and 3.58 x 10" 3 (right). 



routine explores a variety of minima. For small systems 
we obtain identical ground states for all choices of r] con- 
figurations. We conclude that the |r/| = ground state 
represents a robust minimum. For large systems we take 
i] = where convergence is linear in N . We find a variety 
of vortex lattice configurations and mixtures of Mott and 
superfluid- vortex states depending on parameters. In the 
following, however, we focus on slow rotation. 

We now examine the ground state properties of a ro- 
tating system with parameters tuned to give a Mott in- 
sulator at the trap center with a superfluid strip 7 
sites wide) at the edge. For slow rotation (frequencies 
below the Mott gap) the Mott state rotates with the 
lattice giving zero current in the rotating frame while 
the superfluid has non-zero current. The main panel of 
Fig. (fT]) plots the expectation value of the static kinetic 
energy per particle in the rotating frame as a function 
of the rotation frequency. The total static lattice kinetic 
energy measures the net change of the ground state phase 
and therefore drops in steps as the superfluid increases 
vorticity starting from V = 0. The circulation of the 
edge superfluid jumps when the number of effective flux 
quanta passing through the central Mott insulator in- 
creases by an integer to give critical frequencies fiy such 
that: m v /E R w 4V/ (tt 2 |tv| 2 ) . The steps in Fig. (TTJ) 
are slightly parabolic because the hopping term in H 
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FIG. 2: Main panels: The diamagnetic current, J®, (solid 
line) and the edge approximation, J®, (dashed line) plot- 
ted in the rotating frame versus rotation frequency. Insets: 
The superfluid order parameter plotted in the xy plane. The 
top panel is entirely in the superfluid phase (t/U = 0.06 and 
fj,/U = 0.4) while the bottom panel is the Mott phase with 
superfluid edges (t/U = 0.03 and fi/U = 0.4). 
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varies as VI 2 : Re [t exp(iAij)] — t(l 
changes the area of the Mott insulator with fl. 

The change in superfluid circulation can be seen in the 
momentum distribution function and may therefore be 
observable in TOF. The momentum distribution peaks 
associated with the onset of superfluidity expand step- 
wise into rings of radius k r « V / '(|rv|a). The insets of 
Fig. (fT]) show a grey scale plot of pu in the k x k y plane 
for two rotation frequencies. Here we see that a slight in- 
crease of frequency causes a drastic change in the shape 
of the momentum distribution function signaling a jump 
in circulation of the edge superfluid. Observation of the 
number of jumps (i.e. V), k r , and fiy can be used to 
experimentally overdetermine \r^\. 

We now ask if the momentum distribution can yield 
quantitative information related to the edge response. 
The superfluid rotation in the rotating frame can be 
thought of as a diamagnetic current. The top panel in 
Fig. @ plots the excess diamagnetic current, (J D )n — 
(J D )n=o as a function of rotational angular frequency 
deep in the superfluid regime of the trapped BH model 
with t/U = 0.06. The inset shows a grey-scale plot of 
the superfluid order parameter as a function of lattice 
position for Hfl/En = 4.02 x 10 -3 . From the plot we 
see that there is no Mott insulator in the system but 
there is a vortex at the center. The solid and dashed 
lines indicate expectation values of J D and J D , respec- 
tively. In defining the latter we rewrite the parameter 
tv | in terms of an observable, The step corre- 

sponds to the formation of a vortex. Here we see that 
the approximation made in defining J D does not hold for 
bulk current, i.e. (jj 3 ) 56 (Jq > ). The superfluid order pa- 
rameter varies appreciably along the direction transverse 
to the current and, as a result, the diamagnetic current 
cannot be written in the form J . The bottom panel 
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FIG. 3: Schematic of the expected primary momentum distri- 
bution peaks as a function of wavevector arising from clock- 
wise and counter-clockwise propagating edge superfluids sur- 
rounding the Kitaev optical lattice. The diamagnetic contri- 
bution to edge flow splits the peaks in an external magnetic 
field in the Abelian (left) and non-Abelian (right) states. 



shows the same but for a different hopping, t/U = 0.03, 
allowing for a bulk Mott insulator surrounded by edge 
superfluid (see inset). The dashed line reproduces the 
solid line indicating that (J<P) e ~ (Je')e is in fact a 
good approximation for an edge superfluid. Here we find 
a small spatial variance in the superfluid order param- 
eter along the direction transverse to the current, i.e. 
EJ^| 2 (W - \^\) 2 /EM\ 2 < !• We have demon- 
strated, by a realistic simulation, that one can generate 
systems obeying this small variance condition and that, 
as a result, the edge diamagnetic current can be writ- 
ten in terms of an observable, the momentum distribu- 
tion. We propose that in general (J<P) can be restored 
from observation of pk and k and input parameters to 
yield a powerful tool for studying insulating optical lat- 
tice phases with edge states. We now discuss potential 
implications. 

Certain insulators are characterized by a Chern num- 
ber which can be related to their one dimensional edge 
currents 2^, 23}- As an example we assume that the Ki- 
taev model [l7l. Hit [lit can be realized with two compo- 



nent bosons in a honeycomb optical lattice. In Ref. [17 1 
it was shown that the non-Abelian state can be stabi- 
lized in a uniform external magnetic field and that edge 
states exhibit a quantized Righi-Leduc effect. This pre- 
diction asserts that the net edge energy-current displays 
a thermal version of the quantum Hall effect where the 
transverse temperature difference, T, between the bulk 
and exterior of the sample establishes a quantized en- 
ergy current along the edge: jf — vg{irkBT) 2 /6h. Here 
g = 1 for bosons and 5 = 1/2 for fermions. In general we 
expect a clockwise and counter-clockwise energy current 
with an excess number of modes v = v + g — v_q. 

We speculate that, in principle, TOF measurements of 
the momentum distribution function can be used to iden- 
tify chiral edge currents of constituent bosons around two 
dimensional insulators. Flow along the edge of the hon- 
eycomb lattice (+9 and —9) corresponds to concentric 
rings in TOF. If B* is chosen to lie at, for example, the 
V = 1 — > 2 crossing point the +9 and —9 modes will oc- 
cupy different momentum channels resulting in two con- 



centric rings of differing radii in the momentum distribu- 
tion (Fig. $3§). This suggests that, in principle, TOF can 
be used to study chiral edge current and possibly iden- 
tify insulators with non-zero Chern number. In practice, 
however, an observation of edge current in TOF pushes 
current experimental capabilities even for the simplest 
case of a BH Mott insulator. 

Sufficiently accurate observations of wave vector and 
momentum distribution can be used as a quantitative 
probe but are difficult to achieve. Slow rotation induces 
only small k ~ k r modulation of the momentum distri- 
bution peaks. Small features in the momentum distribu- 
tion peaks may not be resolvable experimentally because 
TOF measurements are ultimately limited in k-space res- 
olution [§]. Furthermore, n s can be adversely affected by 
interactions during TOF. Most importantly, the number 
of particles in edge states needs to be sufficiently large to 
overcome background noise in detection. 

We thank CW. Zhang for helpful discussion and ARO- 
DTO, ARO-LPS, and LPS-NSA for support. 
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